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Pade-improvement of four-loop /3-functions in massive </> 4 scalar field theory is shown to predict the known five-loop contribution with 
astonishing (0.2%) accuracy, supporting the applicability of Pade-summations for approximating all-orders MS QCD /3-functions, as 
suggested by Ellis, Karliner, and Samuel. Surprisingly, the most general set of [2|2] approximants consistent with known two-, three-, 
and four-loop contributions to the QCD /3-function with up to six flavours fail to exhibit any zeros that could be interpreted as 
positive infrared fixed points, regardless of the unknown five-loop term. When they occur, positive zeros of such [2|2] approximants 
are preceded by singularities, leading to a double-valued /3-function that is decoupled entirely from the infrared region, similar to 
the /3-function of SUSY gluodynamics. 



Higher order terms of the QCD MS /3-function 
o dx 
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x = a s (fj,)/ir are, upon truncation, known to permit 
the occurrence of fixed points other than the ultravi- 
olet fixed point at x — - 0; e.g. the positive infrared 
fixed point (IRFP) which occurs for 9 < rif < 16 
when the series for fax) in (1) is truncated after two 
terms \fi Q = (11 - 2n//3)/4; fa = (102 - 38n//3)/16; 
xirfp = — fa/ fa]- However, the fixed points arising 
from such truncation are likely to be spurious, as the 
candidate-value for xmpp is sufficiently large for the 
highest-order term in the series faxmFp) to be compa- 
rable in magnitude to lower terms [e.g.|/3iic 3 | = |/3oie 2 |]. 
In a recent paper, u Ellis, Karliner and Samuel predicted 
the coefficient fa via Pade approximant methods, and 
claimed that fa-2 and their prediction for fa yield a 
Pade summation of the /3-function with a nonzero IRFP 
consistent with an earlier prediction by Mattingly and 
StevensonH This Mattingly-Stevenson scenario leads to 
the freezing-out of the coupling to a constant value in 
the infrared region, as shown schematically in Fig. [l]. 

Pade summation of the MS /3-function identifies the 
infinte series 
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{Ri = fa/ fa) with a Pade approximant which incorpo- 
rates the known UV asymptotics of the /3-function, 
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Figure 1: Mattingly-Stevenson scenario 

whose Maclaurin expansion reproduces the known terms 
in the infinite series (2): 
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An IRFP of the /3-function would, in this approximation, 
necessarily be identified with a positive zero of I3\n\m}] 
i.e. a positive zero (x num ) of 1 + a\x + ... + omx n , the 
numerator of Sw\m]> provided SWim] remains positive 
for < x < x num . This latter requirement precludes the 
existence of a positive zero (xden) of the denominator 
1 + bix + ... + 1)mx m that falls in the interval < x < 

Xnum • 

One cannot automatically dismiss the possibility of 
such a denominator zero occuring within the true QCD 
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/3-function. The /3-function of SU(N) SUSY gluodynam- 
ics is known exactly if no matter fields are present ja and 
exhibits precisely such a zero: 



0{x) = 
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1 -Nx/2 
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The /3-function (4) has been discussed further by Ko- 
gan and ShifmanQ If (4) is incorporated into (la), the 
resulting Kogan-Shifman scenario (Fig. |^) for x(fi) is in- 
dicative of both a strong phase in the ultraviolet region 
(the upper branch of Fig. 2) as well as the existence of an 
infrared cut off (/i c ) on the domain of x(/i) that renders 
the infrared region /x < /i c inaccessible. a 
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Figure 2: The Kogan-Shifman scenario [xj en = x(fi c )] 

Without use of Pade summation methods, as de- 
scribed above, the known terms of the infinite series rep- 
resentation (2) for the MS QCD /3-function offer little 
information as to whether the Mattingly-Stevenson (Fig. 
1) or Kogan-Shifman (Fig. 2) scenario is appropriate for 
the evolution of the strong coupling. However, there is 
reason to believe that Pade summation representations of 
the /3-function [eq. (3a)] may indeed be appropriate for 
quantum field theoretical calculations. Ellis, Gardi, Kar- 
liner, and Samueltru have argued that Pade-summations 
(3) converge to their perturbative series (2) as N and M 
increase for any such series dominated by a finite set of 
renormalon poles, consistent with the following asymp- 
totic error formula for the difference between Rn + m+i 
in (2) and the value i?^ 1 ^-^. predicted via use of the 
[N\M] approximant in (3b): I'll 
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a We are assuming a s to be real. The possibility of a s being com- 
plex for fi < fi c is addressed in ref. 4. 
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In (5), {a, b, A} are constants to be determined. 

To demonstrate the utility of this asymptotic error 
formula, consider ±he known /3-function for massive 4 
scalar-field theoryQ 



c= l -(d^){d»4>)+ l -w?4> 2 + g 
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13(g) = l.5g 2 |l - —g + 10.8499 5 2 - 90.5353 ff 3 
+ R i9 4 + ...}. 
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Using the first two terms of (7) to generate a [0 1 1] approx- 
imant, as in (3), one would predict i?^ ade = (-17/9) 2 , 
in which case we see from (5) that 



62 = 



(-17/9) 2 - 10.8499 
10.8499 
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1 + (a + b) 
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Using the first three terms of (7) to generate 
a [1 1 1 ] approximant, one would predict i?4 adc = 
(10.8499) 2 /(-17/9), in which case 



5, = 



[(10.8499) 2 /(-17/9)] - (-90.5353) 
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Equations (8) and (9) are two equations for the unknown 
constants A and (a + 6), with solutions 

A = [1/S 2 - l/Ss}- 1 , (0 + 6) = (S 2 - 2S 3 )/(S 3 -8 2 ). (10) 

We can substitute (10) into (5) to determine R4. 
The first three terms in the series (7) generate a [2|1] 
approximant whose Maclaurin expansion (3b) predicts 
^Padc _ (_90.5353) 2 /(10.8499). Upon substitution of 
A, (a + b), and i?4 adc into (5) we find that this asymp- 
totic error formula predicts that 



i?4 



< ad 7(i + s 4 ) 



(ii) 



i?r ad6 [3 + (a + 6)]/[3 + (a + b)-A] = 947.8 



The value of R4 in (7) has been explicitly calculated^ to 
be 949.5, in very close agreement with (11). The series 
of steps leading from (5) to (11), a methodologicaLrecipe 
first presented in ref. 6, has also been applied Eroto N- 
component scalar field theory, for which the Lagrangian 
(6) is modified such that <f> —> </> a , (j> 2 — > J2a=i ^ > 
N = {2,3,4}. Agreement with calculated values of R4 
(i? 4 = /3 4 //3 ) remains within 3.5% for N < At 

This startling agreement suggests that Pade method- 
ology may also be applicable to the QCD /3-function, 
particularly in the n/ = gluodynamic limit where 
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1. such methods are expected to be most accurate 



2. comparison with the Kogan-Shifman scenario for 
SUSY gluodynmamics is most relevant. 

For rif = 0, the 4-lppp MS QCD /3-function, as defined 
by (la), is given by El 



(3{x) = -^x 2 [l + 2.31818 



41. 53833^ 
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k=4 



The coefficients Rk are presently not known for k > 4. 
The first three terms in the series (12) are sufficient in 
themselves to determine the Pade approximants Su 121 
and Spii], as defined in (3). These approximants are 
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1 - 2.7996a; - 3.7475a; 
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1 - 5.1178a; 



1 - 5.9672a; 



1 - 8.2854a; + 11. 091a; 2 



(13) 



(14) 



In both (13) and (14), the (first) positive denomina- 
tor zero precedes the positive numerator zero: for (13), 
Xnum = 0.264 > Xden = 0.195; for (14), x nam = 0.168 > 
Xden = 0.151. Consequently, x num cannot be identified 
with the Mattingly-Stevenson IRFP in either case, as 
this zero is separated from the small x-region by a sin- 
gularity past which the /3-function switches sign. Indeed 
the ordering < Xden < x num is suggestive of a Kogan- 
Shifman scenario in which x num , if taken seriously is an 
ultraviolet fixed point (UVFP) characterizing the strong 
phase [i.e., the upper branch of Fig 2). 

We can apply the asymptotic error formula (5) to 
the series (12) in precisely the same way we applied it 
to (7). We then obtain an estimate i? 4 = 302.2, analo- 
gous to (11). Using this value of R4 in conjunction with 
the known terms of (12), it is possible to obtain a [2|2]- 
approximant /3-function 
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1 - 9.6296a; + 4.3327a; 2 



1 - 11.9477a; + 23.913a; 2 



(15) 



The first positive numerator zero x num — 0.1092 is 
again larger than the first positive denominator zero 
Xden = 0.1063, precluding the identification of x nU m as 
the IRFP of the Mattingly-Stevenson scenario (Fig. 1). 
Instead, the /3-function (15) is consistent with the Kogan- 
Shifman scenario of Fig. 2, with x num again identified as 
a nonzero UVFP for the strong phase. 

Curiously, the ordering < Xden < x num character- 
izes [2|2]-approximant (3- functions even if i? 4 is allowed 
to be arbitrary. The most general such /3-function that 



reproduces the first four terms of (12) [the first three 
being known] is 
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P[2\2](x) = —£ X X 



(16) 



1 + (13.403 - 0.076215i? 4 )a; - (22.915 - 0.090166iLOa; 2 



1 + (11.084 - 0.076215i? 4 )a; - (56.727 - 0.26685i? 4 )a; 2 

It is easy to verify the first positive numerator zero of 
(16) is always larger than the first positive denominator 
zero [Fig. [|, although these zeros become asymptotically 
close as R4 — > +00. Thus, we see that the first positive 
zero of any [2|2] Pade approximant whose Maclaurin ex- 
pansion reproduces the known terms of eq. (12) cannot 
be identified as an IRFP, nor is such Pade-summation 
indicative of a Fig. 1 scenario for the MS n/ = /3- 
function. 

Remarkably, the same set of conclusions can be 
drawn for the physically interesting case of three light 
flavours. When nf = 3, the 4-loop MS QCD /3-function 
is given byB 
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(3(x) = — — [1 + (16/9)a; + 4.471065a; 2 



20.99027a; 3 + R k ~ 



k=4 



(17) 



with Rk not presently known for k > 4. The known 
terms in (17) determine [2|1] and [1|2] Pade-summation 
representations of the nf = 3 /3-function, 
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1 - 2.91691a; - 3.87504a; 2 



1 - 4.69468a; 



1 



.17337a; 



1 - 9.95115a; + 13.2199a; 2 



(18) 



(19) 



The positive zero of fimi\{x) (x — 0.2559) occurs after 
the pole at x = 0.2130; the positive zero of /3n\i](x) at 
x = 0.1223 similarly occurs after a pole at x — 0.1194. 
The most general [2|2] approximant consistent with (17) 
is 



/3[2|2](a 



9x 2 



(20) 



1 + (7.1945 - 0.10261i? 4 )a; - (11.329 - 0.075643i? 4 )a; 2 



1 + (5.4168 - 0.10261i? 4 )a; - (25.430 - 0.25806i? 4 )x 2 



The Maclaurin expansion of (20) reproduces the series 
in (17), including its (unknown) i? 4 a; 4 term. As was the 
case in (16), the first positive zero of the denominator 
of (20) is always seen to precede the first positive zero 
of the numerator, regardless of R4. Thus the [1|2], [2|1] 
and most general possible [2 1 2] -approximant representa- 
tions of the nf = 3 MS /3-function uphold the ordering 
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< Xden < x num , an ordering that precludes the identifi- 
cation of x num with the IRFP of the Mattingly-Stevenson 
scenario. Moreover, [2|2]-approximant (3- functions for ar- 
bitrary i?4 have been constructed El analogous to (16) 
and (20) for nf = {4,5,6}, and for each of these, the 

< Xden < Xnum ordering persists regardless of R4. A 
range for R4 for which an ordering compatible with Fig. 

1 (0 < x n um < x den) is possible does not occur until 
n f = 7. 

As noted above, the ordering < Xd en < x n um sug- 
gests the occurrence of a double-valued QCD coupling- 
constant, as is the case in SUSY gluodynamics (Fig. 2). 
Such a scenario is seen to decouple the infrared region 
fi < He from the domain of a s , provided a s is under- 
stood to be real. Such a scenario is also indicative of 
a strong phase at short distances u with possible impli- 
cations for dynamical electroweak symmetry breaking, 
suggesting that QCD may even furnish its own "techni- 
colour." 
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Figure 3: Dependence of the first positive numerator zero (x num ) 
and denominator zero (xd en ) of (16) on R4, the (presently un- 
known) five-loop coefficient of the rif = /3-function. 
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